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1 

METHOD OF CONSTRAINING A HETEROGENEOUS PERMEABILITY 
FIELD REPRESENTING AN UNDERGROUND RESERVOIR BY DYNAMIC DATA 

BACKGROUND OF THE INVENTION 

Field Of The Invention 

[0001] The present invention relates to a method of constraining stochastic 
models representing heterogeneous underground zones such as oil reservoirs to 
data referred to as dynamic because they vary with the fluid displacements. These 
data are, for example, production data or pressure data obtained from well tests. 

BACKGROUND OF THE INVENTION Description Of The Prior Art 

[0002] The use of stochastic models of Gaussian type for representing the 
heterogeneity of underground structures is for example described by : 

Journel, A.G. and Huijbregts, Ch. J. : "Mining geostat i st i cs Geostatistics ", Academic 
Press, 1978, or Chiles, J. P. and Delfiner, P. : "Geostatistics - Modeling Spatial 
Uncertainty", Wiley-lnterscience Publishers, John Wiley & Sons, 1999. 

[0003] A numerical reservoir model can be formed from a set of grid cells to 
which the values of a realization of a stochastic model of Gaussian or related type 
are assigned. These values can be assimilated to porosities or permeabilities. 

[0004] Matching the numerical reservoir model with the dynamic data measured 
in the field can be done in form of an optimization problem. A previously defined 
objective or cost function quantifies the difference between the dynamic data 
measured for the real medium and the corresponding responses of the numerical 
reservoir model. These responses are calculated by means of a numerical flow 
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simulator. The goal of the optimization problem is to modify the reservoir model or 
rather the associated realization to minimize the objective function. This process is 
iterative : each iteration implies direct simulation of the flows. A good optimization 
method should allow4e-: a) fRe4tf¥ -modifying realizations discretized on a very large 
number of grid cells; b) carrying out the modifications while respecting the stochastic 
model, jrer that is the modified realization has to be coherent with the stochastic 
model; and c) limiting the number of direct flow simulations because they require a 
considerable calculating time. 

[0005] Simulated annealing can be mentioned as an example of known 
optimization techniques. This approach is for example described by : 

Gupta, A. D. et al. : "Detailed charactorization of fractur e d l i moston o format i on using 
stoch a stic i nv e rs e approach e s Characterization of Fractured Limestone Formation 
Using Stochastic Inverse Approaches ", SPE Ninth Symposium, 1994. 

[0006] This technique is based on realization values e xchang e s exchange 
between grid cells. Upon each exchange, the objective function has to be calculated 
and therefore a direct flow simulation has to be carried out. This process requires an 
excessive number of iterations. Furthermore, in order to preserve the agreement 
between the realization and the stochastic model, an additional term concerning the 
variogram is introduced in the objective function, which makes optimization more 
delicate. 

[0007] Other optimization techniques, more commonly applied, are based on 
gradients calculation. Several approaches based on gradients are presented by : 
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Tarantola, A. : "Inverse probl e m th e ory M e thods for data fitting and mod ol 

param e t e r estimation Problem Theory - Methods for Data Fitting and Model 
Parameter Estimation ", Elsevier Science Publishers, 1987. 

[0008] They require calculating the gradients of the objective function with 
respect to the parameters of the problem which are the values of the realization at 
each grid cell. The realizations are then modified as a function of these gradients so 
that the objective function decreases. The problem related to conditioning of a 
reservoir model to production data is not linear : the minimization techniques using 
gradient calculation are used iteratively. After each modification of the realization, a 
direct flow simulation is carried out and the gradients are recalculated. Applied 
suddenly, the gradient methods lead to calibration of the dynamic data but they 
destroy the coherence between the stochastic model and the realization. Besides, 
they — the gradient methods do not allow to consider a very large number of 
parameters. In order to overcome these limits, geostatistical parameterization 
techniques can be integrated thereto. The pilot point method can be mentioned at 
this stage, which is described by : 

de Marsily, G. et al. : "Interpretation of i nt e rf e r e nc e t e sts in a w e ll field using 
g e ostat i stica l — tochniquos — to — fit the — permeabi li ty — distribution — m — a — res e rvo i r 
mede Unterference Tests in a Well Field Using Geostatistical Technigues to Fit the 
Permeability Distribution in a Reservoir Model " in Verly, G. et al. (ed.), Geostatistics 
for natura l r e sourc e s Natural Resources characterization, Part 2, D. Reidel Pub. Co, 
1984. 
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[0009] This method cons i sts i n s ele ct i ng selects in the realization a certain 
number of points referred to as pilot points, in calculating the derivatives of the 
objective function with respect to the values at these points, in modifying the values 
of these points accordingly and in propagating the disturbance thus defined by 
means of a kriging technique. The pilot point method can induce deviant value 
variations of the pilot points. 

[0010] Another geostatistical parameterization technique, which allows the 
above-mentioned difficulty to be overcome, is the method of gradual deformation of 
a stochastic model of a heterogeneous medium such as an underground zone. It is 
described and used by Hu, L.-Y. et al., in French patents FR-2,780,798 and FR- 
2,795,841 filed by the app l icant assignee . 

[0011] The gradual deformation method allows to gradual l y modif y qradual 
modification of a realization of a stochastic model of Gaussian or related type while 
respecting this model. The deformed realization still is a realization of the stochastic 
model. When the gradual deformation method is introduced in an optimization 
process, the procedure is as follows. The initial realization is combined with a fixed 
number of independent realizations related to the same stochastic model. These 
realizations are called complementary realizations. Combination is controlled by as 
many deformation parameters as there are complementary realizations. It produces 
a new realization. The derivatives of the objective function with respect to the 
deformation parameters are then calculated. The tettei ^deformation parameters are 
modified so as to take into account the information from the derivatives. A first 
optimization in relation to the deformation parameters provides a realization verifying 
the stochastic model and reducing the objective function. In general, this 
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optimization process has to be repeated several times with different complementary 
realizations so as to sufficiently reduce the objective function, which may require, in 
some cases, a prohibitive number of direct flow simulations. 

[0012] The gradual deformation method allows, in some cases, te 
moatf ymodifying a realization locally. This possibility is justified when the gradual 
deformation is combined with the FFTMA geostatistical generator described by : 

Le Ravalec, M. et al. : The FFT moving av o rag o Moving Average (FFT-MA) 
g o norator : an eff i c i ent numerical method for g e nerat i ng and condit i on i ng Gaussian 
simulat i ons Generator: An Efficient Numerical Method for Generating and 
Conditioning Gaussian Simulations , Math. Geol., 32(6), 2000. 

[0013] This generator produces realizations for a stochastic model of Gaussian 
type specified beforehand by convoluting a Gaussian white noise with an operator 
depending on the covariance function. A local deformation can be carried out by 
applying the gradual deformation method to the Gaussian white noise underlying the 
realization. 

[0014] The gradient techniques developed to date for calibration in relation to 
dynamic data are based on a direct link between the variations to which the 
realization representing the reservoir is subjected and the variation of the objective 
function. The optimization process involves modifying first the realization, then 
starting a flow simulation to apprehend the resulting variations for calibration. 

[0015] A different approach valid in cases where the flows are modelled by 
streamlines, is proposed by : 
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Wang, Y. and Kovscek, A. R. : A str e am li ne approach for h i story matching 
production — dat aStreamline Approach for History-Matching Production Data . 
SPE/DOE IOR, 2000. 

[0016] Inverting By inverting the procedure of conventional approaches, tbese-the 
authors propose a methodology based on the successive application of two stages. 
The first stage focuses on calibration : i t a l lows to e va l u a t e allowing evaluation of the 
modification to be applied to the effective permeabilities of the streamlines in order 
to improve calibration. The second stage relates to the transfer of the effective 
permeability variation of the streamlines to the realization. This process is continued 
until a satisfactory calibration is obtained. The main limit of this approach is that it 
does not preserve the coherence of the realization in relation to the stochastic 
model. 

SUMMARY OF THE INVENTION 

[0017] The method according to the invention allows to work working out 
numerical models, representing heterogeneous underground media such as oil 
reservoirs or aquifers, in accordance with a set of dynamic data measured in 
production, injection or observation wells, and with a stochastic model characterizing 
the spatial distribution of the heterogeneities. 

[0018] tt- The method essentially comprises minimization of the objective function 
from an iterative process subdivided in two stages and by integration of the gradual 
deformation method in the second stage. The first stage identifies the disturbance to 
be applied to the effective permeabilities of predetermined zones so as to reduce the 
objective function. The second stage focuses on the minimization of the difference 
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between the desired effective permeabilities and the corresponding effective 
permeabilities calculated for the realization considered, by modifying this realization 
by means of the gradual deformation method. This minimization stage is an 
intermediate optimization which requires no new flow simulation and allows te 
defofm deforminq the realization while respecting the Gaussian stochastic model. If 
the objective function is not satisfactory at this stage, the process is resumed at the 
first stage. 

[0019] The method according to the invention is a decoupled gradual optimization 
method. The term gradual applies insofar as deformation of the numerical model is 
carried out according to the fundamental principles of the gradual deformation 
method. The term decoupled r e m i nds means that the optimization is based on two 
stages, the first er>e ~staqe defining the correction to be brought to the effective 
permeabilities, the second stage propagating l^ -the correction independently to the 
numerical model. These two stages imply that the flow simulations are carried out 
only during the first stage in order to improve calibration of the effective 
permeabilities. On no account do e s t The second stage does not require a flow 
simulation to estimate the way the numerical reservoir model should be deformed. 

[0020] In other words, the method according to the invention allows te 
fefm forming iteratively a model representative of the permeability field of a 
heterogeneous medium, discretized by a grid, such as an underground zone, 
constrained by a priori geologic data and dynamic data collected in the medium by 
means of measurements and observations obtained beforehand. I^ -The method 
includes : 
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a first stage comprising : 

a) generating an initial model of the permeability in accordance with a 
Gaussian or related stochastic model, coherent with the a priori geologic data, and 
carrying out, by means of a suitable simulator, a simulation of the fluid flows, 

b) identifying zones inside the reservoir, calculating the effective 
permeabilities of these zones and, from the simulator results, estimating the 
corrections to be brought to these effective permeabilities to improve calibration in 
relation to the data, and 

a second stage comprising : c) propagating the corrections to the whole of the 
grid cells of the permeability field, by means of an iterative optimization process 
comprising minimizing a functional which quantifies the difference between the 
effective permeabilities required to obtain said— the calibration and the effective 
permeabilities calculated for the permeability field considered, using a technique of 
gradual deformation of realizations of the stochastic model. 

[0021] The zones are defined either manually or automatically from the flow 
simulator. 

[0022] According to an embodiment, simulation of the flows is carried out by 
means of a streamline simulator and the zones of the medium are identified by a set 
of grid cells traversed by one or more streamlines of fixed geometry. 

[0023] According to another embodiment, said- the zones are identified as volume 
portions on the periphery of wells running through said- the medium, within the 
framework of well tests. 
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[0024] According to another embodiment, at least one gradual deformation 
parameter is assigned to each zone. 

[0025] In relation to the prior art, the method according to the invention allows te 
ca li brat e calibration of a realization of a stochastic model to dynamic data while 
keeping the coherence in relation to the stochastic model and by reducing 
significantly the number of direct flow simulations to be carried out. Unlike the 
approaches developed to date, a deformation of the realization does not 
systematically involve a new flow simulation. J^ The method thus affords the 
possibility of better exploring the space of the realizations and of rapidly determining 
not only one constrained realization, but several ones. Furthermore, the method 
according to the invention allows to d o form deforminq the realization from a large 
number of parameters, which accelerates the intermediate optimization. It also 
allows deformation of the realization by zones, and the latter zones can be the same 
as the zones predetermined for calculation of the effective permeabilities during the 
first stage. 

Brief Description Of The F I GURES Drawinqs 

[0026] Other features and advantages of the method according to the invention 
will be clear from reading the description hereafter, given by way of non limitative 
application example, with reference to the accompanying drawings wherein : 

Figur e Fig. 1 shows an example of an optimization scheme in two stages 
according to the method of the invention , and 

F i gur e Fig. 2 shows a comparison between a real fractional flow and the 
corresponding simulated fractional flow. 
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Detailed Description Of The Preferred Embodiments 

[0027] Traditionally, calibration of reservoir models is an iterative process 
wherein a flow simulation is carried out each time the reservoir is disturbed. For 
simplicity's sake, we cons i der that the numerical reservoir model is considered to 
com e s come down to a realization of a stochastic model for the permeability. In 
other words, a flow simulation is required for any variation of the permeability field. 
On the contrary, the method according to the invention can propose several 
permeability field variations using a single flow simulation. J VThe method according 
to the invention therefore consists of has two stages (Fig._1). During the first stage, 
zones of the reservoir are determined and the correction to be brought to the 
effective permeabilities of these zones in order to improve data calibration is 
estimated. Then, an optimization process is started to propagate the disturbance 
determined for the effective permeabilities of these zones to the permeability field 
representing the reservoir. The first stage requires a flow simulation b e caus e it 
which depends on a comparison between real data and the corresponding synthetic 
responses. On the other hand, the second stage requires no additional flow 
simulation. Finally, the optimization performed in the second stage is parameterized 
according to the gradual deformation method, which allows to pr e serve preserving 
the coherence of the permeability values distribution (or any other property 
considered) in relation to the spatial variability model. 

[0028] The general algorithm is summarized as follows : 

An initial permeability field in accordance with a stochastic model being 
generated, 
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a) stage 1 - a flow simulation is performed, the effective permeabilities of 
these zones are calculated, zones are identified in the field which is considered and 
the variation to be applied to the effective permeabilities of these zones in order to 
reduce the difference between the production data measured in the field and the 
simulated production data is estimated; 

b) stage 2 - the disturbance required at the level of the effective 
permeabilities of the zones is propagated to the whole of the permeability field by 
means of an optimization process. The function to be minimized quantifies the 
difference between the desired effective permeabilities and the effective 
permeabilities calculated for the permeability field considered. The gradual 
deformation technique is used to modify the permeability field; 

c) return to b) as long as the calibration is not satisfactory. 

[0029] The effective permeabilities of the zones can be calculated using a 
simulator or thev- the effective permeabilities can be obtained by means of averaging 
techniques well-known to th e man skil le d in the art. 

[0030] The method according to the invention is adaptable to any flow simulator 
insofar as means to define effective permeabilities for different zones of the reservoir 
are available, the latter being identified manually by the user or automatically from a 
given criterion. For example, for a well test simulator, rings of increasing radius, 
c e ntr e d centered on the wells, can be selected to define the zones. The effective 
permeabilities can in this case be related to the apparent permeabilities. 
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[0031] In the example developed hereafter, w e focus on the case of a streamline 
simulato r is the focal point . The zones considered here are determined by the 
streamlines themselves. Each one of the aforementioned two stages is described in 
detail in the two sections hereunder. 

Stage 1 Modification of the e ff e ct i v e — p e rm e ab i l i t ie s Effective 
Permeabilities 

[0032] The objective of the first part is to describe the geometry of the 
streamlines and to identify the modification to be applied to the effective permeability 
of these lines so as to better calibrate the production data. A line or a set of 
streamlines define a zone. 

[0033] Modelling flows by streamlines involves certain characteristics. The 
displacement of a fluid along a streamline is a one-dimensional problem; the 
streamlines do not communicate with one another. When the mobility ratio of the 
fluids i n pr e senc e is one and when the boundary conditions do not change, the 
geometry of the streamlines is fixed. When the mobility ratio is different from one, 
there are two a lt e rnat i v e s alternatives : the geometry of the streamlines is fixed and 
the flow variations during the fluid displacement is allowed or the flow is uniformly 
distributed between the current tubes streamlines and the geometry of the 
streamlines is periodically updated. We — choos e — tThe first configuration 
straightaway is chosen . 

[0034] To simplify the problem, w e consid e r groups of streamlines are considered 
rather than the streamlines ind i v i dua ll y individually : we4Thus r e duc e J he number of 
parameters is reduced . The effective permeability of a group of streamlines is 
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expressed as the harmonic mean weighted by the fluid volumes of the permeabilities 
of the grid cells traversed by the streamlines : 

Kf = Eq.l 

k=\ 7=1 kj 

[0035] Kj eff is the effective permeability of the i-th group of streamlines. N g si_ is the 
number of streamlines contained in this i-th group, N k is the number of grid cells 
traversed by the k-th streamline of the i-th group, qig is the flow for the k-th 
streamline at the level of the j-th grid cell. Atkj is the flight time for the k-th streamline 
through the j-th grid cell. 

[0036] For a given initial realization, the flow simulation allows to d e term i n o 
determination of the geometry of the streamlines and to calculate calculation of their 
effective permeabilities. The effective permeability variation that would bring a 
production data calibration improvement remains to be evaluated. The fractional 
flows observed and the corresponding simulated fractional flows (Fig._2) are 
therefore compared for the producing wells. The streamlines are then arranged in 
the increasing order of their breakthrough times and the fractional flow curves are 
discretized. The segments Aq thus defined are associated with groups of 
streamlines. In the absence of accordance between the simulated flows and the real 
flows, the effective permeabilities of the groups of streamlines are considered to be 
responsible for the differences. A correction applied to all of the streamline groups 
allows this difference to be reduced. Consider the group of streamlines associated 
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with a flow increment Aqj (Fig._2). The effective permeability desired for streamline 
group i in order to improve calibration is : 



T^eff i* simulated jy-eff 

^i, desired ~ ^ i, simulated 

i, reference 



where preference and ^simulated are the mean breakthrough times associated with 
increment Aqj for the real and simulated fractional flow curves. The same procedure 
is carried out for each group of streamlines. 

Stage 2 : Intermediate opt i m i zation Optimization - Modification of the 
perme a bi li ty Permeability fiel dField 

[0037] The first stage determines a variation to be applied to the effective 
permeabilities of the streamlines to calibrate the fractional flows. The goal of the 
second stage is to transfer this variation of the streamlines to the permeability field 
while preserving the coherence of the permeability field in relation to the stochastic 
model. An optimization problem is therefore defined. The function referred to as 
intermediate function is to be minimized here : 



1 N 8 

r 1 V^i, desired ^ i simulated ) 



[0038] N g is the number of streamline groups. Kf desired is the desired effective 

permeability for the group i of streamlines in order to reduce the differences between 
the measured and simulated fractional flows : this value has been determined during 
the first stage. Kf simulated is the effective permeability of the group i of streamlines for 
the permeability field considered. This optimization problem is non-linear and can 
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involve as many parameters as the permeability field contains grid cells. By 
integrating the deformation method as the parameterization technique, w e can 
r e duc e t he number of parameters can be reduced and provide a permeability field 
modification that matches the stochastic model defining the spatial distribution of the 
heterogeneities in the permeability field. Furthermore, in order to be able to modify 
the permeability field by zones, the zones corresponding to the various streamline 
groups, the gradual deformation is not applied to the permeability field itself, but to 
the underlying Gaussian white noise (Fig. 1 ). 

[0039] The simplest version of the gradual deformation method cons i sts i n 
comb i n i ng combines two Gaussian white noises, Zi and z 2 , according to the relation : 

z(x)(p) = z x (x)cos(np) + z 2 (x)sin(;rp) 

where p is the deformation parameter and x the position. The gradient of z with 
respect to p is : 

— = -m x svnfjtp) + nz 2 cos(7tp). 
dp 

[0040] The Gaussian white noise z is then converted to a Gaussian realization y 
of covariance model C, average zero and variance one, by means of a convolution 
product developed at the heart of the FFTMA generator : 

y = f*z. 

[0041] f results from the decomposition of the covariance function. The derivative 
of y with respect to the deformation parameter is : 



V 
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dp dp 

[0042] Any other geostatistical generator producing Gaussian or related 
realizations can be used in place of the FFTMA generator, provided that calculation 
of the gradients is integrated therein. The FFTMA generator however affords the 
advantage of fast execution, even for realizations discretized on a very large number 
of grid cells. 

[0043] Realization y, which is for the moment c e ntr e d centered and reduced, 
f^ vthat is of average zero and variance one, can then be converted to a Gaussian 
realization w of average m and variance a 2 : 

w(x) - m + oy(x). 

[0044] The derivative with respect to the deformation parameter becomes : 

dw dy 
—(x) = a-f-(x). 
op dp 

[0045] At this stage, static data observed at precise points, in wells for 
e xamp le s example . also have to be taken into account. This information is in general 
integrated in the realization generated by means of a kriging technique. The 
constrained realization w c is deduced from : 



™ c (*) = W dK (*) + 0(*) - W K (*)) 
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where w d K and w K are the realizations obtained, for the first one, from kriging of the 
real observations and, for the second, from kriging of the values of w at the 
observation points. The kriging estimator, in the dual frame, is expressed as follows : 

n 

1=1 

[0046] C is the covariance function. The Xj are the positions of the n observations. 
The weights pi are independent of the position, but they depend on the deformation 
parameters. W e It can be shown that the derivative of w with respect to the 
deformation parameter is obtained from : 



dw„ , . dw dw K 



dp dp dp 



dw. 



dp ttdp 
% dp dp 



[0047] The permeability field k is deduced from the lognormal transformation of 

w c : 

k(x) = exp(w c (x)). 

[0048] The permeability gradient with respect to the deformation parameter is 
expressed as follows : 



dp dp 
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[0049] Equation 1) allows 



^ calculation of the effective permeability of a 



group of streamlines for a given permeability field. The derivative of the effective 
permeability with respect to the deformation parameter is deduced therefrom : 



[0050] These various relations show how to deform a permeability field and to 
calculate the gradients of the effective permeabilities of the streamlines with respect 
to a deformation parameter. All these relations can be readily generalized in cases 
where several deformation parameters are involved. One can notably decide to 
assign a deformation parameter to each zone defined for calculation of the effective 
permeabilities. 

[0051] Several techniques can be considered for minimizing the intermediate 
function, but since the gradients of the effective permeabilities are available for an 
insignificant calculation cost, a Gauss-Newton type approach may be judicious. To 
determine the disturbance to be applied to the deformation parameters so as to 
reduce the intermediate function, the following system is solved : 

HAp = -V(F/). 

[0052] Ap is the disturbance to be defined, V(FI) comprises the gradients of the 
intermediate function with respect to the deformation parameters and H is an 
approached matrix of the Hessian matrix : 



dKf 



dp 



eff k=l J= l 




k) dp 
N*.; Ar *.; 



H = G WG. 
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[0053] G is the sensitivity matrix : tt -which includes the derivatives of the effective 
permeabilities of the streamline groups with respect to the deformation parameters. 
W is the weight matrix 4-fc- which is here equal to the identity matrix. 

[0054] Finally, the algorithm developed to minimize the intermediate function is 
described as follows. Stage 1 allows to d e f i n e defininq , on the one hand, the desired 
effective permeabilities for the streamline groups and, on the other hand, zones 
assigned to these groups. 

a) At least one deformation parameter is assigned by zone. 

b) An initial Gaussian white noise and at least one complementary 
Gaussian white noise are defined. 

c) A gradual deformation is carried out by combining the previous 
Gaussian white noises. The derivatives of the Gaussian white noise z resulting from 
the gradual deformation with respect to the deformation parameters are 
simultaneously calculated. 

d) The Gaussian white noise z is converted to a permeability field k and 
the derivatives of the permeabilities with respect to the deformation parameters are 
also determined. 

e) The effective permeabilities of the streamline groups defined in stage 1 
and their derivatives are calculated. 

f) The disturbance to be applied to the deformation parameters to reduce 
the intermediate function is then estimated. 
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g) At this stage, several alternatives emerge. If the intermediate function 
is not weak enough and if it does not appear to have converged, the deformation 
parameters are updated and ewe-returns to c) occurs . If the intermediate function is 
not weak enough but seems to have converged, ene-returns to b)_pccurs, irer that is 
the initial Gaussian white noise is updated and a new complementary Gaussian 
white noise is randomly selected. Finally, if the objective function is sufficiently weak 
or if it appears that it is not going to decrease further, stage 2 comes to an end. 

[0055] The procedure described here focuses on the case of streamline flow 
simulations. Streamlines actually appear to be a very natural tool for defining zones. 
The ~From the h e art of the present invention involving identification of zones and 
calculation of the effective permeabilities for these zones, this choice seems to be 
logical. Other types of application can however be considered. For example, the 
information relative to the various zones could be provided by well tests. The 
effective permeability can in this case be calculated as a function of the radius of 
investigation around the well ; the zones considered are then rings. The flow 
simulator used for well test simulation can be a standard flow simulator. The 
procedure to be followed is similar to the procedure described for the streamlines. 

Flow simulation allows to idefrtifv identification of zones and te 

determ i n e determination of the effective permeability for these zones, which can be 
compared with the data measured in the field. Then, minimization of an intermediate 
objective function according to the approach described above allows to propagat e 
propagation the correction to be applied to the effective permeabilities to the 
absolute permeabilities of the grid cells in the zones while respecting the a priori 
spatial variability model. 



21 

ABSTRACT 

A Method method of constraining a stochastic model of Gaussian or related type, 
representing a porous medium such as an underground reservoir, to data 
characteristic of the displacement of the fluids is disclosed, r 

The method is based on an iterative development in two stages. The first ene 
consists i n carryinq stage carries out a flow simulation, in identifying zones inside the 
reservoir and in estimating the modification to be brought to the effective 
permeabilities of these zones so as to improve calibration between the real data and 
the corresponding responses obtained with the flow simulator. The second stage 
involves an intermediate optimization problem intended to minimize an intermediate 
objective function (implemented from a gradual deformation technique) measuring 
the difference between the effective permeabilities calculated for the zones and the 
effective permeabilities identified during the first stage for better calibration. 



